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ABSTRACT. We analyze the potential for liner buckling in a 40-in Kevlar49/epoxy 
overwrapped spherical pressure vessel (COPV) due to long, local depressions or ‘valleys’ in 
the titanium liner, which appeared after proof testing (autofrettage). We begin by 
presenting the geometric characteristics of approximately 20 mil (0.02 in.) deep depressions 
measured by laser profilometry in several vessels. While such depths were more typical, 
depths of more than 40 mils (0.02 in.) were seen near the equator in one particular vessel. 

Such depressions are largely the result of overlap of the edges of overwrap bands (with 
rectangular cross-section prepreg tows) from the first or second wrap patterns particularly 
where they start and end. We then discuss the physical mechanisms of formation of the 
depressions during the autofrettage process in terms of uneven void compaction in the 
overwrap around the tow overlap lines and the resulting 10-fold increase in through- 
thickness stiffness of the overwrap. We consider the effects of liner plastic yielding 
mechanisms in the liner on residual bending moments and interface pressures with the 
overwrap both at the peak proof pressure (-6500 psi) and when reducing the pressure to 0 
psi. During depressurization the Bauschinger phenomenon becomes very important 
whereby extensive yielding in tension reduces the magnitude of the yield threshold in 
compression by 30 to 40%, compared to the virgin annealed state of the liner titanium. In 
the absence of a depression, the liner is elastically stable in compression even at liner- 
overwrap interface pressures nominally 6 times the - 1000 psi interface pressure that exists 
at 0 psi. Using a model based on a plate-on-an-elastic-foundation, we develop an extensive 
analysis of the possible destabilizing effects of a frozen-in valley. The analysis treats the 
modifying effects of the residual bending moments and interface pressures remaining after 
the proof hold as well as the Bauschinger effect on the compressive yield threshold. The key 
result is that depression depths of up to 40 mils can be tolerated, but above 40 mils, the 
Bauschinger effect drives destabilization, and buckling becomes increasingly likely 
depending on the details of depression formation during autofrettage. It is almost certain 
that destabilization and buckling will occur for depression depths beyond 55 mils. The main 
equations and formulas for treating the various phases of depression development and 
potential buckling, are only briefly outlined in the paper, but are available from the authors. 

I. Introduction 

P roof testing of composite-overwrapped pressure vessels (COPVs) that have a metal liner typically result in 
plastic yielding of the liner material so that upon returning to zero pressure, the liner is in a state of hoop 
compression under an overwrap-liner interface pressure. The main benefit of this autofrettage effect is that 
during COPV service under steady pressure, which may involve cycling from zero pressure to maximum operating 
pressure, the liner stays within the elastic state and fatigue life is greatly extended. At the same time, placing the 
liner in a compressive state when unpressurized introduces the risk of local buckling. Even if unbuckled after the 
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proof test, the potential for liner buckling may remain in marginal vessels should an unpressurized vessel with a 
carbon or Kevlar/epoxy fiber overwrap be raised in temperature, since the overwrap and liner have opposite 
coefficients of thermal expansion. Although designing against liner buckling is ostensibly straightforward, several 
factors may increase the risk such as variability in liner and overwrap mechanical properties, bonding between the 
liner and overwrap (if required for stability), and irregularities in tow and band packing. 

In this paper we develop an analysis of the potential for liner buckling in a 40-in diameter, Kevlar49/epoxy 
COPV with a 0.74 inch thick overwrap and 0.105 inch thick titanium liner. Some vessels were found to have local 
depressions or valleys (dents) in the titanium liner that appeared after proof testing. We begin by presenting the 
geometric characteristics of approximately 20 mils (0.02 in.) deep depressions measured on the inside liner surface 
by laser profilometry (and seen as bumps). While such depths are more typical, depressions of more than 40 mils 
(0.02 in.) were seen near the equator in one particular vessel. These depressions are believed to be the result of 
overlap of the edges of the overwrap bands from the first wrap patterns particularly where they begin and terminate. 
The physical mechanisms of formation of the depressions during the autofrettage process involve uneven void 
compaction and overwrap stacking at the locations where bands overlap, especially at the beginning and ends of 
wrap patterns. We also discuss the consequences of extensive liner plastic yielding mechanisms on residual bending 
moments and interface pressures both at the peak proof pressure of ~ 6500 psi, and during depressurization to 0 psi 
after autofrettage. In the latter case the Bauschinger phenomenon on lowering the yield threshold becomes very 
important since the extensive yielding in tension reduces the magnitude of the yield threshold in compression by 30 
to 40%, compared to the virgin annealed state. 

To provide a context for the detailed analysis and modeling, we performed a preliminary analysis of the in-plane 
buckling loads and critical half-wavelength scales from various results from classical buckling theory, including (i) a 
free-standing sphere under external pressure, (ii) a short cylinder under a uniaxial compressive load, and (iii) a plate 
on an elastic foundation with hinged boundary conditions. By comparing the results of these various cases we found 
that the analysis of the potential for buckling of the liner in these liner/overwrap configurations reduces to 
considering the behavior of a rectangular hinged plate (a 0.105 in. thick liner section spanning zero bending moment 
lines) on an elastic foundation (provided by the 0.74 in. thick overwrap) and with a “frozen-in” depression or valley 
that also has residual bending moment and interface pressure fluctuations remaining after plastic yielding during 
autofrettage. In the absence of a depression, we find that the liner is elastically stable in compression at liner- 
overwrap interface pressures nominally six times the approximately 1000 psi that exists at zero vessel pressure after 
autofrettage. The analysis treats the modifying effects of the residual bending moments and interface pressures 
remaining after the proof hold at peak pressure, as well as the Bauschinger attenuation effect on the plastic yield 
threshold in compression. The main equations and formulas for treating the various phases of the depression 
development, and potential buckling during depressurization, are only briefly outlined in the paper, but are available 
from the authors in more detail. 

First we discuss the fundamental differential equations and their solution for transverse liner elastic deformation 
due to application of a narrow perpendicular line load. This is done both for the case of in-plane tension, as occurs 
during autofrettage, and in-plane compression, as occurs during subsequent depressurization. Using independently 
determined material stiffness parameters, we calculate depression or valley profiles that match well those 
experimentally observed using laser profilometry. A crucial aspect is that the transverse (through-thickness) 
stiffness of the overwrap in this initial pressurization phase is governed by overwrap porosity from extensive voids, 
and thus, is an order of magnitude less that the transverse stiffness that exists after void collapse and during 
depressurization. This key aspect is the reason why permanent volume growth from autofrettage varies by more 
than a factor of two among these vessels. 

Next, we determine mathematical results for interface pressure and residual bending moment profiles that remain 
after liner yielding is complete at the maximum proof pressure, which is held for a few minutes. An important 
aspect is the calculation of the extent of dispersion or spreading of the sharp spike in interface pressure, which 
results from liner yielding and flattening of the overwrap bump. Also important is the fact that force balance 
requires that the integral of the pressure fluctuation must sum to zero so that the pressure peak must be compensated 
by pressure reductions on each side of the pressure hump of the displacement valley. Through two integrations 
across the locally parabolic pressure peak, the dispersion smoothes out the pointed profile of the bending moment (a 
characteristic of line loads such as in three-point bending), replacing it with a locally parabolic shape. This 
dispersion has virtually no effect on the depression profile shape since this shape involves two further integrations of 
the bending moment in the underlying fourth order differential equation. However it smoothes out the pressure 
fluctuations and make it easier to appreciate the local interface pressure loss with increasing depression depth, which 
is a precursor to buckling. Within physically reasonable limits, the actual breadth of the dispersion has very little 
effect on the main predictions regarding the onset of buckling. 
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Third, we calculate the changes in the interface pressure and bending moment profile from the state at peak proof 
pressure, as well as changes in the depression shape and depth that occur during depressurization of the COPV from 
the peak proof pressure down to zero pressure. The key driving shape in the equations is the depression shape for 
the tensioned liner during the pressure buildup of autofrettage, which occurs at a low transverse stiffness of the 
overwrap. However, in the depressurization, transverse overwrap stiffness becomes the high stiffness of the 
compacted overwrap with voids removed. Finally we discuss the key framework for superimposing the results for 
various stages of the autofrettage cycle, particularly the depressurization phase. Also considered are criteria that 
allow final assessment of the likelihood of liner buckling for various values of initial ‘frozen-in” depression depth. 

The key result is that depression depths of up to 40 mils can be tolerated by the liner. Above 40 mils, the 
Bauschinger effect dominates and destabilization and buckling becomes increasingly likely depending on the details 
of depression formation during the proof pressurization. Whatever the mechanical property variations for a vessel, it 
is almost certain that destabilization and buckling will occur for depression depths beyond 55 mils. 

II. Laser Profilometry Measurements and Potential Buckling Mechanisms 

From laser profilometry on the inside surfaces of both 26 and 40 inch Kevlar/epoxy COPVs at zero pressure, such 
valleys (referenced from the overwrap side), have been seen in certain vessels. These depressions typically have a 
characteristic width of a few inches along the liner surface and in some cases tend to exhibit a modest overshoot on 
each side, i.e., small hills flanking the valley, suggestive of mild oscillatory behavior in the disturbance as it dies out. 
An illustration of such a depressions found in a 26 in. spherical vessel is shown in Figure 1, where the horizontal 
axis is circumference measured in degrees and the vertical axis is degrees latitude relative to the equator at 0 degrees 
[1]. The smaller ripples have width 1.26 inches, consistent with the width of the bands in the wrap patterns. 

circumferential degrees around equator 
200 250 300 350 

I | | I 


Figure 1: Profilometry map of surface depression (blue) observed in a 26 in. Kevlar/epoxy COPV with a 
titanium liner. The depression, measured at NASA-JSC WSTF [1] was 22 mils deep. 

Figure 2 shows a similar depression observed in a 40 inch spherical Kevlar 49/epoxy COPV. A somewhat 
surprising feature is that the ripples have spacings much more consistent with the bands in the 2 nd , and higher wrap 
patterns than in the 1 st wrap pattern. The former are about 1.65 inches wide while the latter are 1.05 inches wide. 
An interesting issue is whether the small ripple wavelength is more consistent with and somewhat excited by a 
natural buckling wavelengths of the overwrap-liner combination. Nonetheless the span of the main depression 
appears to be significantly larger than of the periodic ripples. 

Figure 3 shows the dimensional characteristics of a depression. The liner is put into compression when the vessel 
pressure is reduced to zero after autofrettage, and thus, a critical question arises regarding the possibility of 
buckling, and its dependence on the depth and width of the valley, as well as other characteristics of the liner and 
overwrap. Since the liner is not bonded to the overwrap (in fact the surface has a release agent by design), the 
possibility exists that the liner will buckle if the depression depth formed at proof, exceeds some critical value 
whereby the liner loses the elastic support of the overwrap and separates. Of primary interest are the conditions 
under which such liner buckling can occur, and “safe” conditions under which the liner will be stable. 
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Figure 2: Profllometry map of a surface depression (red region) observed in a 40 in. diameter Kevlar/epoxy 
COPY with a titanium liner. The depression, measured at NASA-JSC WSTF [1], was 20 mils deep. 



Distance along liner, inches 

Figure 3. Characteristics of the cross-section of an approximately 20-mil depression in a vessel [1]. 

A. Mechanisms of Initial Formation of the Liner Depression 

In modeling the effects of such a liner depression, the first issue to consider is the conditions under which a 
depression forms as well as the “frozen-in” geometric bending moment and interface pressure characteristics it 
inherits from autofrettage. Immediately after curing the freshly wound vessel and before autofrettage, there are no 
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obvious mechanical reasons why such depressions should occur. However during the pressure build-up of the 
autofrettage cycle, which includes a 2 to 5 minute hold at 3000 to 3500 psi and a 2 minute hold at the proof pressure 
of 6520 psi, the conditions are favorable to form such depressions. Interpreting experiments on through-thickness 
deformation of vessel cross-sections it is found that, the transverse compressive overwrap stiffness through the 
thickness before autofrettage is up to a factor of 10 lower than it is after autofrettage (e.g., about 25,000 psi versus 
250,000 psi). The apparent cause of this large increase is permanent collapse of air voids in the overwrap as the 
interface pressure increases during autofrettage. 

From plots of growth in boss-to-boss length versus pressure, there is strong evidence that much of the void 
collapse already occurs before or at the hold pressure of 3000 to 3500 psi during autofrettage, even before the tensile 
yield strain in the liner has been reached. An important observation is that the permanent volume growth resulting 
from the proof process is from 1.5 to 2.25 times larger that calculated from both FEA and membrane mechanics 
models of the vessel proof process, again the likely result of compression of entrained air voids. This discrepancy 
translates into liner radial displacements of about 40 mils beyond predictions, or about 6% of the 0.74 in. overwrap 
thickness. An uneven void distribution in the overwrap is likely where overlap of bands in the first wrap occurs as 
they travel from the equator to the boss zone, particularly in starting and terminating wrap patterns and circuits. 
This makes it possible for local bumps in the overwrap to occur as voids on each side are collapsed. Coincidentally 
depressions of 20 mils, as shown above in Figures 1 to 3, are approximately the thickness of two stacked, 
Kevlar/ epoxy prepreg tows made from 4560 denier yarn. Stacking of multiple overlaps is possible where wrap 
patterns begin and end if special effort is not made to ensure offset positioning. Thus depressions deeper than 20 
mils are possible, and a depression of 40 mils was observed in one vessel near the equator. 

An important aspect is the permanent cross-sectional shape of such “frozen-in” depressions or valleys in the liner, 
including depth, width and length. Extensive liner yielding in tension occurs during proof and the final plastic strain 
at the end of the 2 minute proof hold at 6520 psi is larger than the elastic strain at yield onset, by 15 to 40% (i.e., 
total strain is 2.15 to 2.4 times the yield strain). This relieves most of elastic bending moment induced by the 
change in liner curvature as the depression forms. The residual bending moment and associated residual interface 
pressure disturbance is critical to assessing the likelihood of buckling as the vessel pressure is lowered to zero and 
the liner is put in compression. This compression results in approximately a 1000 psi average interface pressure (in 
both 40 inch and 26 inch diameter COPV)) between the overwrap and liner, which is crucial to the liner’s stability as 
the overwrap acts as an elastic foundation. If this interface pressure locally falls to zero for larger depression depths, 
elastic foundation support is lost and the conditions become favorable for liner buckling. 

B. Potential Mechanisms of Buckling and Inherent Resistance Factors 

Without the elastic foundation support of the overwrap, the liner is not capable of supporting the approximately 
1000 psi interface pressure that exists after autofrettage when the vessel returns zero pressure. As we show later, the 
critical buckling pressure of the stand-alone liner is at most 600 psi with a critical half wavelength of about 2.5 
inches. Thus to prevent buckling, the liner relies on the lateral support of the overwrap as an elastic foundation. 
Maintaining a positive interface pressure is crucial to this support mechanism, since the liner and overwrap will then 
act as though they are bonded. An appreciation of the effectiveness of the overwrap in preventing liner buckling can 
be obtained from analysis of a plate on an elastic foundation, as we develop shortly. We find that the reinforcing 
effect is by a factor of about 6 meaning that the liner could nominally sustain interface pressures of close to 5,000 
psi (assuming it did not yield in hoop compression) as part of the coupled liner-overwrap system, and the half 
wavelengths of buckling would be about 1 inch. The fact that the liner is not bonded to the interface means that the 
likelihood of liner buckling is closely tied to the likelihood of loss of interface pressure and subsequent liner 
separation from the unbonded overwrap. By itself, this local loss of interface pressure may not be sufficient for 
buckling. It is also necessary that characteristic length over which contact has been lost, and the distance between 
zero bending moment points, must also match the characteristic buckling length of an unsupported plate-like region 
of the liner with hinged boundary conditions. From buckling theory for free-standing plates, it can be shown that for 
in-plane compressive loads that exist at zero vessel pressure (about 10,000 lb./in), the characteristic, half-wavelength 
for buckling is about 2 inches. Thus a model for buckling of the liner must consider the likelihood such conditions 
will be realized when the depression depths and widths are of sufficient magnitude. 

As mentioned, a crucial aspect of when such conditions might arise is the magnitudes of local residual bending 
moments and interface pressure increases that result from the formation of such a depression during autofrettage. 
As the depression forms there is extensive liner yielding that occurs that will relax the natural bending moments 
from liner curvature changes, whether elastic or plastic. It can be shown that for a laboratory section of the liner in 
three-point bending the maximum elastic bending moment possible (before yielding on the surface begins at about 
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128 ksi) is about 240 lb-in./in. and the ultimate plastic bending moment is 3/2 times this value, or 360 lb-in. per 
inch. However these large magnitudes are only possible if the in-plane liner loads average to zero, which is not the 
case during autofrettage where the in-plane strain of the liner becomes more than double the elastic yield strain 
when the depression formation process is active. Thus more than half of the elastic bending moment from the 
axisymmetric, triangle-shaped bending stress distribution is flattened and eliminated (all on one side and part of the 
other side), and for deeper depressions, the reduction is expected to be about 70%, since the bending is likely to 
begin earlier. This stress relaxation occurs even for small induced moments, which are approximately proportional 
to the depression depth for a given width. 

C. The Bauschinger effect on Attenuation of Compressive Yield Threshold 

As the vessel pressure is lowered from the peak proof pressure to 0 psi pressure, further bending will occur that 
changes the curvatures relative to those that are already frozen-in as a result of the pressure build-up and yielding of 
proof. This bending will superimpose new triangular stress distributions through the liner thickness on top of the 
residual blunted ones remaining from the peak proof pressures. Depending on the location in the depression, the un- 
blunted portion that was at a lower tensile stress than the tensile yield stress 128 ksi, may become severely blunted 
as the vessel pressure is reduced to zero and the liner is put into compression. This blunting occurs even though the 
nominal or average in-plane compressive stress generated upon unloading to zero pressure is only -80 ksi. The true 
cause of this premature blunting phenomenon is the Bauschinger effect, whereby the effective yield stress of the 
liner material when placed into compression, after having suffered a large yield strain in tension, becomes 
approximately -80 ksi and not -128 ksi as was found in tests performed at NASA Glenn Research Center [1]. 

Uniaxial load cycles applied to a specimen taken from a Ti liner 



Strain (in/in) 

Figure 4. Bauschinger effect in attenuating the compressive yield stress after autofrettage 

The most obvious consequence of the Bauschinger effect is a substantially reduced bending moment compared 
the elastic bending moment associated with such curvature changes from the frozen-in state. The magnitude of the 
bending moment reduction depends on the details of the bending history and yielding during the void compaction of 
autofrettage as well as the location in the depression where the bending moment is evaluated. The magnitude of the 
liner bending moment in the valley area that can safely be sustained at zero vessel pressure is not likely to be much 
more than 30 percent of the maximum elastic bending moment observed in a simple three-point bending experiment, 
or about 72 lb-in./in. During depressurization of some vessels, even if the paths of bending moment versus liner in- 
plane compressive load allow for significantly higher bending moment limits at zero vessel pressure, the 
Bauschinger effect leads to softening in liner bending response reflected in a moderately reduced secant modulus 
This will have profound consequences on the potential for liner buckling. We estimate that this will lead to a 
collapse in interface pressure and thus liner separation from the overwrap. This condition quickly becomes the most 
critical condition for buckling of the liner. 
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III. Preliminary Buckling Analysis 

To study the formation and effects of a liner depression or valley, initially in terms of liner tension during 
autofrettage and subsequently in terms of liner compression upon vessel depressurization, we consider an analysis of 
a plate on an elastic foundation potentially modified to account for spherical curvature effects. We also must 
consider cylindrical curvature since the depressions are approximately long ellipses in projected outline onto the 
surface of the sphere as shown in Figures 1 and 2. Thus they are long enough for ring stiffening effects to play some 
role. We first consider some classic results to frame the problem. 


A. Classical Buckling of a Sphere and Simple Correspondence to Plate Buckling 

We let E, be the liner tensile modulus , v, its Poisson’s ratio, t, be its thickness and R : be its radius. In the 
classical theory for buckling of a sphere, the critical external pressure, P cnl , and associated critical in-plane (hoop) 
load , N cril , and buckling half- wavelength, JL cnt /2 , are respectively 


P = 

cnt 


2 E,tf 




P R 

_ crit l 


E,tf 


\'/ 4 


R 




Kn,l 2 = ^^l 




(i) 


For the liner of the 40-inch Kcvlar/cpo xy COPV we have t, =0.106 in., R t =19.19 in. , E t =16.5xl0 6 psi and 
v ; = 0.342 so for the liner itsMf under pressure we obtain P cnt = 619 psi , N crit = - 5,936 psi and X cnt /2 = 2.50 in. 
Thus for the OMS liner alone with no overwrap support, this result predicts collapse of all liners at zero vessel 
pressure, since the interface pressure is about 960 psi, which is well above the calculated 619 psi. 

Later we develop the theory for a plate (liner region) on an elastic foundation (the overwrap), but certain results 
for a simply supported plate are mentioned here. In Euler buckling, the critical, in-plane buckling load, N crit plau , , for 
a simply supported square plate of length and width, L, and under biaxial loading is 


N, 


cr it, plate 


x 2 E,t] 
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Taking the spherical buckling half- wavelength for L, i.e., L = X crit / 2 from (1), yields 
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Surprisingly, the flat-plate plate buckling load is the same as that for the in-plane load of the sphere when its half 
wavelength is the same as the plate dimension, L. This is a key observation, because we can study more 
complicated buckling of a composite sphere by replacement of critical regions with flat plate sections. In particular 
we can study the effect of the elastic foundation support of the overwrap on the behavior of the spherical liner. 


B. Key Equations for Buckling of a Plate on an Elastic Foundation 

The fundamental differential equation for the transverse displacement v(x) of a plate on an elastic foundation in 
terms of axial distance, x , is the 4 th order equation 


E,t i d 4 v(x ) 

12 ( 1 - 1 x 7 ) dx 4 


d 2 v[x ) 
dx 2 


+ — v(x) = 0, 
to 


-00 < x < 00 


( 4 ) 


where E oi is the overwrap transverse modulus (through-thickness compressive stiffness), t 0 is the overwrap 
thickness, N t is the "running” in-plane liner load (with a negative value in compression) and the remaining 
parameters are as defined in connection with (1). This differential equation for transverse (indentation) 
displacement, v (x), can be written in terms of a non-dimensional load, N, and stiffness parameter, <D , as 


d 4 v d 2 v 

— ^-2N — r- + <Dv = 0, — 00 < x < 00 

dx 4 dx 2 


N = 


6N,(l-v ? ) 
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E/t, t 0 


( 5 ) 
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C. Possibility of Reinforcement Effects due to Spherical Shape compared to Flat Plate 

Before proceeding, we digress briefly to consider the relative importance of curvature effects associated with the 
surface contours of the liner depressions, which are roughly elliptical in projected outline on the sphere surface. 
One might expect the radial curvature along the major axis of the elliptical indentation (approximately the curvature 
of the liner itself) to act to suppress the simply supported plate buckling along the minor axis (narrow direction of 
the depression), which will suffer major curvature changes associated with protruding into the inside of the vessel. 
The theory for cylinder buckling developed by Sam Batdorf can be used to shed light on the importance of this 
effect. Regarding Batdorf s theory for a short cylindrical shell (with no overwrap support), and with hinged end- 
conditions, one can derive his results beginning with the liner-only differential equation 

e£ 


d 4 v d 2 v E,t, 
-tf»-rr + -£-v = 0 


12^1 -v, 2 ) dx 4 1 dx 2 R 2 
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and adding to this the elastic foundation support effect, which results in 
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In the spirit of obtaining the dimensionless equation (5) we define 
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For the second stiffness term in parentheses to provide a contribution comparable in magnitude to the first, one 
would need E 0t « E, t,t 0 j R 2 , and for the 40-inch COPV under study this requires E o t « 3510 psi , which is much 
too small (by almost two orders of magnitude as we see later). On the other hand, the natural buckling load, N , , 
(with only the cylindrical support) and the natural half-wavelength, /L d / 2 , are found, respectively, to be 
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= 2.50 in. 


( 10 ) 


Surprisingly, these are exactly the same results as (2) and (3) for the sphere above. Thus when considering buckling 
of the liner depressions or indentations, we need only consider plate-like effects in the depressed regions and assess 
the extent to which overwrap support will suppress buckling or its absence will allow buckling to occur. 

We also calculate the buckling half-wavelength for a flat, square, plate-like region (or cylindrical ring) with zero 
bending moments at each end and free edges on the side (no edges in the ring), and subjected to the in -plane load 
N, p=0 = -9,172 lb/in. , the value at 0 psi in the 40 inch COPV. Later we show this is the in-plane liner load at 0 psi 
vessel pressure. This half-wavelength happens to be the classic Euler result, which for the 40-inch vessel gives 


Kff/2 


7T 2 E£ 


|6(i-wR^ c 


= 2.00 in. 


( 11 ) 


This establishes a critical length-scale of 2.0 inches between zero bending moment points acting as hinges. If such 
lengths exceed 2.0 inches, buckling becomes a risk particularly if the overwrap support becomes ineffective once the 
interface pressure falls to zero, indicating the onset of liner separation from the overwrap. 


IV. More Detailed Buckling Analysis 

Next we consider a more detailed analysis, considering both the pressurization stage during autofrettage when 
overwrap compaction and through-thickness stiffness increase occurs, and liner yielding and indentation occurs 
from tow and band cross-over pile-ups. Then we consider the depressurization phase where further liner yielding 
and depressurization may occur if the indentation is sufficiently deep. First we present the solution to the main 
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differential equation (5) subject to the boundary conditions associated with an imposed displacement, S b , or 
equivalent load along a line perpendicular the in-plane load, N, and an imposed zero slope along that line. (The 
case of non-zero slope has also been solved but is not presented.) 

The solution begins with the usual complex exponential form for such equations, namely, v = Ce sx , which 
yields the characteristic equation, s 4 -21% 2 +<D = 0 , which has the squared roots 


/V<[> - N 2 , <D > N 2 
a/n 2 -®, N 2 > <t 


( 12 ) 


A. Liner in In-Plane Tension. 

When N is positive and the liner is in tension, we could have the second case where s 2 is real provided N is 
sufficiently large. However this condition requires N > a/o , and for E 0l = 250,000 psi in the compacted state 
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= 50,100 lb/in 


(13) 


but for the 40-inch vessel, the maximum value for N, in tension is the liner yield value, N, = a, y t, = 13,568 lb/in , 

so s 2 cannot be real. This means that the disturbance due to overwrap protrusion will have a decaying sinusoidal 
character with mild overshoots (rather than pure exponential decay). This is consistent with the profilometry results 
described earlier. Thus we need only consider the case <J> > N 2 and this separated into two sub-cases 


f N, N > 0 (liner in tension) 

1 — |N|, N<0 (liner in compression) 


(14) 


The case of loading in tension is important to the formation of the disturbance shape during autofrettage 
( N > 0 and <J> - N 2 ). In this case characteristic equation ,? 4 - 2Ns' 2 + <T> = 0 factors as 


(s 2 - 4®e +w ) (s 2 - J®e~ w ) = (.? + ^e +w ' 2 ) (s - </®i e +w ' 2 ) (s + V® e~ w/ 2 ) (s + V®e~ w/2 ) = 0 


where tan# = -N 2 )/n 2 , thus yielding four roots, s l9 s 2 , s 3 and s 4 given by 

,s | 2 = a/o (cos (0/2) ± /sin ((9/2)) , s 34 = — a/o (cos((9/2)±;sin((9/2)) 

Since tan 2 9 + 1 = l/cos 2 6 we obtain cos# = 2cos 2 (#/2) -1 = l-2sin 2 {0/2) = 1 - n/a/o so that 


cos 


m)=\ 


l + cos# 


I VO+N 
V 2a/o 


sin(#/2) = 


1-cos# 


/ a/O-N 

V 2a/o 


Thus we have the general solution form 

v(x) = C 1 e-^ cosW2)4 cos(^sin(#/2)x) + C 2 e-^ sinW2)4 
+ c 3 e^ cos{l>l2)x cos(^/® sin (#/2) x) + C ^ sHe ' 2)x 


sin (VOsin(#/2)x) 
sin ^ a/o sin(#/2)xj 


(15) 


(16) 


(17) 


(18) 


which applies separately on each side of the line displacement(load) which we place at x = 0 . There can be no 
unbounded solutions as x — > ±oo so that C 3 = C 4 = 0 for x > 0 , and Cj = C 2 = 0 for x < 0 . 

For the boundary conditions we prescribe a displacement and slope, respectively, as v(0) = -# A and 
dv(0)/dx = 0 . Using absolute value functions in x as well as symmetry and continuity, we obtain the displacement, 
slope and curvature, respectively, as 

v(x) = S b e~'^ m ^ el2 ^ |cos|a/® sin(#/2)xj + cot(#/2) sin( J Vo sin(#/2)|x|j|, -oo<x<oo (19) 
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-oo < x < oo 


( 20 ) 


dv(x) \AF -</®cos(»/2ite' • l t[T- . I Q I~\ \ 

and 

= £ 6 ^/FA^ _v ^‘ :os(<^,/2, '' jcos(/AF sin ( 0/2 )xj- cot ((9/2) sin | ^/o sin(0/2)|x|)J, -oo < x < oo (21) 

The depression shape half wavelength is computed by noting that the frequency of the oscillation is 
/ = -x/cp sin(0/2)/(2;r) and the half-wavelength is thus 


K ! 2 


K 

sin (0/2) 


( 22 ) 


This is the “decay” oscillation wavelength of the depression, and it might be matched to profilometer zero crossing 
points of displacements. However this is misleading since in studying experimental profiles one should really 
consider the locations of the zeros of v(x) in order to determine points relevant to characteristic lengths. In (19) 

we note that for x > 0 , the zero crossing for x > 0 actually occurs at -\AF sin (0/ 2)x + 0/2 = n so that an effective 
half-wavelength (width of the valley) might be measured as double this length, or 


^■d.eff / 2 


271-0 
/AF sin (0/2) 


(23) 


From the second derivative or curvature once can calculate the points, x = ±x m 0 , of zero bending moment. 

From (21), these can be shown to be x m 0 = ±(0/2)/|-V® sin(0/2)j . The key point is that the zero points for the 

bending moment are very different from those of the displacement itself. The total length equivalent to an effective 
characteristic buckling half-wavelength, associated with the first point of zero bending moment on each side, is 
double the above value in (24) or 


4,o/2 


0 

Z/O sin (0/2) 


(24) 


Also, this length is different from the width of the valley itself, (23), or the mathematical wavelength, (22). This 
means that modeling the depression as a half-cosine wave with half-wavelength X d eff j2 or even X dep j2 can lead to 

errors in the bending moment zero points thus missing the possibility of the depression turning into a buckle if 
separation begins to occur at the interface. 

B. Liner in In-Plane Compression at Zero Vessel Pressure 

The solution approach in this case is similar to the case in liner tension just treated. For the boundary condition 
y b = 0, 8 h * 0 , the displacement, slope and curvature are respectively 


v(x) = -S b 
0v(x) 

= O, 

dx 


g -t/®sin(9/2|i 

COS (0/2) 
cos(0/2) 


cos^xAF cos (0/2) |;t| - 0/2 j , -oo < x < oo 
-V®sm(0/2)|i| s j n |4/(j COS (0/2) xj, -00 < X < 00 


and 


dv 1 (x) _ I 


— \-d- = S b - — e — — — — cos(-\AFcos(0/2)x + 0/2), — oo < x < oo 

dx cos (0/2) v ' 

Here we note that the mathematical half wavelength of the depression decay is 


(25) 

(26) 


(27) 


-W 2 


n 

<AF cos (0/2) 


(28) 
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( 29 ) 


The effective half- wavelength based on the zero-crossings of the displacement, v(x) , is 

; h 9 + n 

d ' effl v/o cos(#/2) 

and for the zero bending moment location an effective half-wavelength is 
mo/ cos(d/2) 


(30) 


Again this apparent half-wavelength resulting in buckling concern will be much smaller than the apparent half- 
wavelength of the waveform itself given by (28). 


C. Analysis of Indentation and Deformation during Overwrap Compaction Developed During Autofrettage 

In addition to the frozen-in displacements due to autofrettage, there is a frozen-in bending moment distribution, 
which is not proportional to the second derivative of the frozen in displacement, because plastic yielding of the liner 
occurs during autofrettage that greatly reduces this moment especially at the point of largest curvature at the deepest 
part of the liner depression. The process of generating the frozen-in displacement profile comes mainly from void 
compaction when the liner is in substantial tension and the liner overwrap interface pressure is elevated because of 
the high proof pressure. The relevant tension, N, c , may be less than the peak value during autofrettage, denoted as 

N, . Pressure versus boss-to-boss displacement records and experiments on overwrap sections show that much of 
this compaction occurs well before the liner tension reaches N ; , thus it is appropriate to consider the sensitivity of 
the results to any differences between N , c and N, . 

Regarding the frozen-in, bending moment distribution, since the peak vessel pressure during the proof hold of 
autofrettage (about 6520 psi in a 40 inch vessel) is much higher than the point where yielding begins (beyond which 
N,=N, ), the plastic yield strain at the end of the proof hold is more than double the elastic strain when yielding 

begins (at about 3200 psi to 3400 psi). Thus the bending moments at the peak proof pressure will be less than half 
the values calculated on the basis of bending of a purely elastic liner (ignoring yielding), and if the depression is 
sufficiently deep, the peak bending moment is reduced to about 30% of the nominal elastic value. However, the 
bending moment distribution at this point, where N, - N, , is different from what it becomes at the reference point 

N/ = 0 , passed through during depressurization. Likewise, the associated distribution of liner-overwrap interface 
pressure is also different and some care is necessary to distinguish between them. 

We begin by calculating the bending moment distribution and associated interface pressure distribution at the 
peak proof pressure point where N,=N, and the peak displacement is -8 b . We then follow the unloading 

process from the peak point down to the reference state, N, =0 (at ~ 2650 psi). This will allow us to compare 
displacements and arrive at an appropriate reference displacement, -8 b . Once this has been accomplished, we 
can follow the various deformation components, bending moments and interface pressures all the way from the 
maximum tension at autofrettage N l = N t , down through the point of zero tension, N l = 0 on vessel 

depressurization to the maximum compression at 0 psi vessel pressure, N ] = -N t u . 

After winding and curing, and before autofrettage begins, the liner surface is flat and there is no interface 
pressure at all (except for a small effect of the winding tension). However as the autofrettage process begins and the 
vessel pressure is increased from 0 psi, compaction occurs in the overwrap. The effect on the liner displacement 
relative to an equilibrium state (whereby the liner radius increases on average as it physically expands into the 
overwrap) can be viewed as a local depression in the liner that develops relative to this equilibrium state around a 
line perpendicular to the in-plane tensile load. This depression can be viewed as the disturbance that results from a 
very narrow, but very stiff region of the overwrap with many fewer voids and which undergoes very little 
compaction. The physical cause is likely the result of sudden overlap of the edges of two adjacent bands in the same 
wrap pattern as they proceed from the equator to the boss. In general overlap must occur in virtually all bands of a 
wrap pattern as the circumference, at increasing latitude, decreases relative to that at the equator, such that the bands 
can no longer lie beside each other in parallel. Each band must therefore jump and cross over the band laid down 
just before it. This is the cause of the rippling effect seen in Figures 1 and 2 . However each wrap pattern must have 
a starting point, and the first band laid down will have a different overlap circumstance since it must, in fact, be flat 
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everywhere since there is no previous band for it cross over. Note that at bands are stacked four deep since each 
wrap pattern has two passes around the vessel of 360° each and full coverage of the vessel occurs with just 
1 80° degrees of a pass. The result is that, the total stacking of tows in the bands at the beginning and end of a wrap 
pattern and at some latitude location, may be two higher than on average around the rest of the vessel. This will 
cause a local bump in the overwrap and consequently a depression in the liner. 

To capture the effects of this phenomenon, the non-dimensional differential equation describing the transverse 
liner deformation from an equilibrium position is nominally given by 


dx 4 


--2N 


dx' 


- + cp p v p = 0, -oo < X < 00 , 


N,= 


6N,' P (l-v?) 

E/tf 




i2£ o ,. 0 (i - -yf ' ) 

Efi tp 


(31) 


The relevant transverse stiffness of the overwrap for forming the depression shape is the low value E 0 1 0 and not 
E ot , which applies during depressurization after full overwrap compaction and the transverse stiffness is about 10 
times larger. Also the tension is taken as N \ rather than N, c since the effects of any differences are small. 

As the autofrettage process proceeds with increasing vessel pressure, a local depression develops in the liner 
around a line perpendicular to the in plane tensile load. This depression can be viewed as being caused by a 
displacement imposed along a line caused by a corresponding line force from an infinitesimally narrow, region that 
undergoes very little compaction compared to elsewhere in the overwrap during the pressure buildup of autofrettage. 
The effect is nominally captured through the boundary conditions 


/ , dv(x) 

^(0) = -^, hm 7 =0, 


dx 


lim v (x) = 0 

W— »co P V ' 


(32) 


where S b is the depression depth that finally results at the peak proof pressure when N, = N, p . 

The solution to (31) and (32), determined earlier in slightly different notation, can be written as 

v p (x) = -8 ip e^ x \ {cos(fti 0 x)+^ 0 sin(® 0 |x|)}, -oo<x<oo (33) 


where in the present case 

= V®7 sin (A/ 2 )> cos(6> 0 /2) = 


/®, +N , 


2 J® 


> sin(6> 0 /2) = 




2 JO 


<Po = 


+N „ 


®,- N , 


Also several useful identities can be demonstrated, namely 

sin 2 (0 o /2) = ^-, cos 2 (0 o /2) = -^-, ®„=ft> 0 4 (l + % 2 )\ K p =co;{<pl-\) 

cp'+l qj '+ 1 

Thus the solution to the governing differential equation (31) and its four derivatives reduce to 

v p (*) = -5 h p e-‘°°' Po ^ {cos(® 0 v)+^ 0 sin(® 0 |x|)j, -oo < x < oo 

dv ( x ) 

— p - = sgn (x)8 bp co 0 (l + sin(® 0 |x|), -oo < * < oo 


dx 
d\{x)_ 
dx 2 

d\{x) 
dx 2 


= S h p o)g (\ + cpl^e {cos(ffl 0 jc)-p 0 sin(® 0 x )J, -oo < x < oo 


■ - - sgn 


(x)2S bp col<p 0 (l + % )ef‘ wH \ cos(co 0 x) 


2 % 


sin 


(®oW) 


and 


.(*) 


dx 4 


= - 4 S b,M<Po ( 1 + % 2 ) s ( x ) + ^ p O p e~ ahnlxl |cos (ai 0 x) + cp 0 sin (cu 0 |x| 
+ 2S b , p (°l ( 1 + % 2 )®o {<Po -l)e"^ H {cos(ai 0 v)-% sin(® 0 \x\)) 


= -4^®oVo (i + <Po )<?(*) ~ ®o V P (*) + 2N P 


.(*) 


dx' 


(34) 

(35) 

(36) 

(37) 

(38) 

(39) 


(40) 
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where <5(x) is the Dirac delta function. Thus the governing differential equation for studying the depression 
nominal characteristics becomes 


d 4 v„ 


(x) 


-2N 


d 2 v„ 


(x) 


®p v p (*) = - A5 i, P 0J l% {^ + ( Pl)S{x), -oo < x < oo 


(41) 


dx 4 p dx 2 

subject to no other boundary conditions other than v p (x) — » 0, |x| — > oo . As stated this is a nominal solution, but 

modifications are required to obtain a more realistic description that accounts for the dispersive effects of liner 
yielding during the peak tension of autofrettage. 


D. Residual Bending Moments and Spring-Back Pressures at the Liner-Overwrap Interface. 

As the peak pressure of autofrettage is approached, liner yielding takes place that relieves much of the bending 
moment and resulting stress in the liner in the depression region. This means that the pressure distribution will be 
more dispersed than reflected by the concentrated scaled line force, -4 8 b p co\(p a (l + (p \ ) 8 (x) , on the right-hand side 

of (41). Nonetheless a significant residual bending moment will remain and can be viewed as generating a spring- 
back force that is manifest as an interface pressure disturbance relative to the nominally uniform interface pressure 
over the liner-overwrap interface. We assume this liner yielding limits the maximum obtainable bending moment at 
the peak deflection to M f max = K pf max M e max where M e is the maximum elastic bending moment (beyond which 

plastic yielding occurs) one would obtain in a simple three-point bending experiment, and 0 < K pf max < 1 is the 
fraction that would remain if the specimen were subsequently (while maintaining its bent shape) strained in tension 
to some value more than twice (about 2.2 times) the uniaxial yield strain. A reasonable value is K pf max «0.3. 

Beyond some distance away from this peak, say |x| > |x| , we assume the yielding generally lowers the bending 
moment to a value K p times the nominally induced value, where K p » 0.45 . However, within the region close to 
the peak, denoted -|x| < x < |x| , we assume the bending moment is driven by the form of the local but somewhat 
relaxed and diffused pressure distribution around the original line load (or displacement). This local pressure 
distribution, denoted, p p (x) , will be taken to be parabolic (2 nd order polynomial) with 4 th order edge terms to 

provide continuity and a smooth first derivative at the transition points x = ±|x| . This pressure disturbance will 

integrate twice to result in a 6 th order polynomial for the bending moment profile over the region - |x| < x < |x| , but 

will actually be dominated by a 2 nd order parabola rather than the triangular peak in (38) resulting from the line load 
(as naturally occurs in three point bending), which nulls out. The resulting displacement, which requires two more 
integrations, is virtually indistinguishable from the displacement given by (36). However, the details are important 
in determining an accurate description of the pressure distribution, which is critical for determining the possibility of 
liner separation from the overwrap triggering buckling, as the vessel pressure decreases to zero and the liner is 
forced into in-plane (hoop-type) compression. 

We begin with the original form of the equation under the concentrated line load, whereby the overall interface 
pressure profile is obtained from (40) and (35) as 

Pune {*) = F p (*) + K p EI \ p <4 (! + <pl ) 2 {cos (® 0 x) + <p 0 sin (co 0 |x|)} , - oo < x < oo (42) 


where the flexural rigidity, line load, and a convenient line-load like quantity ( force per unit length) are given by 

EI = F p (x) = -K p EI48 Kp a>l<p 0 {\ + (pl)8(x) and F p = K p EI4S bp co } 0 (p 0 (l + cpl ) (43) 

v > / 

Thus the local pressure distribution is 


Pp( x )=~ f p 5 ( x ) + 




-e {cos(ai 0 x)+^ 0 sin(oi 0 x)j, -oo<x<oo 


(44) 


Due to the effects of yield relaxation and compaction, we assume the line load -F 8 (x) becomes diffused and 
can be replaced by the pressure profile described by the 4 th order polynomial 


13 

American Institute of Aeronautics and Astronautics 



Pp( x ) = -P™.[}- A { x lx) 2 +5 (*/*) 4 )> -*<*<* ( 45 ) 

where (-x,x) is the zone over which the p p (x) acts. We assume that both the magnitude of the pressure and the 
slope are zero along the transition lines x = +x , that is, p p (±x) = dp p (x)/dx | + , = 0 . These two equations result in 

Pp( X ) = - Anax (i - 2 i X / X ) 2 + ( X / X ) 4 ) (46) 

Also, we require force balance with respect to the original line load so that J p p (x)dx = —F p from which we get 
/’max = (15/8) F p /( 2x) and thus the diffused local pressure profile is 


Pp ( x ) = 




0 , 


-x < x < x 
— oo < x < — x or x<x<oo 


(47) 


and overall the pressure profile is 

P p (*) = P„ (x) + (F p co 0 ( <p 2 0 +l)/(4 {cos(ft> 0 x)+^ 0 sin(® 0 |x|)}, -oo < x < oo , (48) 

Next we consider the shear profile around x = 0 and first calculate the shear distribution from the local 
distributed pressure p p (x ) over -x<x<x. This is S p (x) = -J p p (tj)dtj which yields 




(49) 


Calculating the total shear profile in the beam requires care. We know that for the line load, the shear profile is 

F, 


Sme (*) = K „ EI ~ = -sgn(x)^e'^ w (cos 


(® 0 x)— ^ — -sin(® 0 |x|H, — oo < x < oo 


(50) 


= F occurs at x = 0 . For the distributed version of the load the overall 


where a jump of magnitude 2^F p j2^ 
decrease of magnitude F occurs more gradually over -x < x < x , which is accomplished by way of 


M*) = 


-sgn(x)-^e 


-sgn(x)^ft 


COS 


(® 0 x)- 


2 <P 0 


sin 




X < \x < 00 


cos 


(® 0 x)--^ — -sin(® 0 |x|)-l| + 5 p (x), 0 < |x| < x 


(51) 


This profile and its derivative are both continuous so the pressure profile is also continuous as was initially assumed. 

To obtain the moment distribution, we note that M p (x) = K p EI d 2 v p (x)/(7x 2 = j S p ( x)dx and adapting the 
results for a line load given earlier, whereby we obtain 

M iine{ x ) = K p EId2v p( x )/ dx2 = (^ J ,/(4^o)) e "‘ W,H {cos(ft> 0 x)-% sin(® 0 |x|)}, — oo < x < oo (52) 


we find that 


M p (x) = 


__p — e | cos ^® o x)-^ 0 sin(® 0 |x|)j, x<|x|<oo 

F 

M p (x) + — — - — |cos(<y 0 x)-^ 0 sin(® 0 |x|)J, 0 < |x| < x 


(53) 


whereby 
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( 54 ) 


M p (x) = | S p (x)dv + sgn(x) — J \dv 


Note that M p (0) reduces to 




/ 

r / \ 2 1 

1 

/ \ 4 1 

1 

i- - 

1- - 

2 

V 

w _ 

6 

w _ 


M,{ 0)= ^ 


4® 0 % V 8 





(55) 


However, we limit M (0) to the maximum value M pf max and this allows us to solve for x using (54)) yielding 


5« 0 % 


x M P f, max 4 ^ 


(56) 


E. Deformation during Depressurization due to the Frozen-in Depression Occurring during Autofrettage. 

Here we develop an analysis of deformation and potential buckling of a plate on an elastic foundation as a result 
of a frozen-in, transient wave disturbance that has occurred autofrettage and is characterized by behavior at the peak 
proof pressure and maximum liner tension. This disturbance is residual deformation induced from a sharp bump 
from fabrication (band crossovers and tow pileups), but modified by overwrap compaction and plastic yielding in 
the liner. Earlier we considered the effects of a residual bending moment and associated disturbance in interface 
pressure also with reference to the peak proof pressure and associated maximum liner tension. On depressurization 
of the vessel the newly induced disturbance from the compressive load in the liner may have a characteristic 
wavelength that is different from the one induced by the original sharp disturbance during pressurization and the 
bending wavelengths may be much longer and different in shape from the frozen in shape. The main cause is the 
greatly increased transverse or through-thickness stiffness of the overwrap from void compaction. In this subsection 
we derive a mathematical description of the behavior on depressurization. 

We let the initial, frozen-in, transverse displacement be of the form 

v 0 (x) = -S bo e~‘°° n ^ jcos(<y 0 x) + ^ 0 sin(cu 0 |x|)j, -oo<x<oo (57) 


where 8 b 0 is the frozen in displacement corresponding to a reference state with zero liner tension, and a> 0 and <p 0 , 
with subscripts ‘0’ given by (34), are waveform parameters associated with the initial, frozen-in, geometry, which 
were appropriate during the pressurization. The displacement parameter here is taken as 8 b 0 since the reference 

liner tension state is different from that for 8 h at the proof peak, though their magnitudes are nearly the same. 

Note that the original model parameters without subscripts, <D , N , were defined in (5) in terms of geometry and 
mechanical properties of the plate and elastic support of the overwrap and are used now to refer to the 
depressurization phase of the vessel since the relevant through-thickness stiffness of the overwrap is E a t , the high 

stiffness of the compacted state. The parameters N and O ; in (31) were used to refer to the tension phase during 
autofrettage where the eventual “frozen-in” profile was generated under the low overwrap through-thickness 
stiffness, E o t 0 , prior to compaction. Thus N ; 0 is in-plane, tensile liner load appropriate to the autofrettage phase 

at peak pressure and A, , a variable, is in-plane liner compressive load upon depressurization after autofrettage. 
Thus the basic model parameters without subscripts are given in (5) and with subscripts are 


N„ = 


6^,o 

E,t] 


and 


®o = 


12(1 ~vf)E 0jS 

E/tj to 


= <t, 


Also 


(58) 


cos(6>o/2) 




tj/o7cos(6> 0 /2) and (p 0 


tan(0 o /2) 


(59) 


and appreciate that E OIO « E 0l after overwrap compaction, the difference being up to an order of magnitude. 
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We let the total transverse displacement, v , be the sum of the incremental displacement, v , induced by the 
scaled applied axial load, N , and the ‘frozen-in’ displacement, v 0 , i.e., v(x) = v(x) + v 0 (x) so that 


v(x) = v(x)-v 0 (x) 


(60) 


is the transverse displacement difference caused by the loading. On depressurization, the differential equation for 
the transverse displacement of a plate section of the liner and overwrap is the 4 th order equation 

^ - 2N ^ + <Dv = 0, — oo < x < oo (61) 

dx 4 dx 2 


Note that the moment of the applied axial load depends on the total displacement (hence no ‘over-tilde’ on the 
second derivative), whereas the transverse displacement and bending components depend on the displacement 
increment. From (60) we can write (61) as 


d 4 v „^ T d 2 v 


— -2N^- + Ov = 2N^-, 
dx dx dx 


-oo < x < oo 


(62) 


where the boundary conditions now are v(x) — > 0 as x — > ±oo and dv{0)/dx = d 3 v{0)/dx 3 = 0 and we note that 
v(x) = 0, -oo<x<oowhen N=0. In terms of the actual displacements, v(x) , the key differential equation is 


d 4 v d 2 v d 4 v 0 

— --2N — - + d)v = -r- + ®v„, — oo < x < oo (63) 

dx 4 dx 2 dx 4 0 


where the boundary conditions become v(x)— >0 as x— >±oo, and dv(0) / dx = d 2 v{0)/ dx 2 =0. However, it is 
simpler to solve the problem in terms of the incremental disturbance upon unloading, v , and then to add the initial 
frozen-in disturbance, v 0 , to get the final result, v(x) = v 0 (x) + v(x) . 

We will need several derivatives of the initial displacement profile (57) to solve the problem. We consider the 
case, x > 0 and from symmetry arguments construct the full solution for -oo < x < oo . Thus from (57) we have 


dv 0 

dx 


S b 0 co 0 {<pl + 1) e ° w sin (oi 0 x) , x > 0 


(64) 


and 


d\ 

dx 2 
d\_ 
dx 3 


= 8 bfi a>l{<pl+\)e °* >nx |cos(ry 0 x)-^ 0 sin(r» 0 x)}, x>0 
= S hfi al{(pl +l)e“ w |(^ 0 2 -l)sin(® 0 x)-2^ 0 cos(® 0 x)}, 


x > 0 


(65) 

(66) 


~^r = ^,,o®o (^o + l)e~°* nx {(^ 0 2 + 2(p 0 - i) co {(o$c)-<p 0 {% -2cp 0 -l)sin(<y 0 x)|, x>0 (67) 

Thus the right-hand side of the key differential equation (62) is 

jl 

2N 2 ° = A b e~ Dx {©, cos(ft> 0 x) + © 2 sin(o 0 x)} (68) 

where 

A i = 2N 5 b 0 wl {(pi + 1) , D = w Q (p a , ©, =1 and 0 2 = -(p a (69) 

Thus for x > 0 , our key equation to solve is 

74 ~ i2 ~ 

— — -2N — ^- + <59 = A b e~ Dx {©, cos(<» 0 x)+© 2 sin(r» 0 x)} (70) 

dx dx 

The four complementary solutions from solving the homogeneous equation (left-hand side set equal to zero) are 
the same in structure as obtained in solving (5), and given by 

Vj = e~°’ ipx cos(ftix) , v 2 = e~ mpx sin (ox) , v 3 = e mpx cos(<yx) and v 4 = e 0 ’ 9 ’ 1 sin(®x) (71) 
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We also require a particular solution, v(x) = v (x) , to (70) and by the method of undetermined coefficients we can 
work with the form 

v p = e“ Dr (,Pcos(<y 0 x) + (2sin((y 0 x)) (72) 

where P and 0 are constants to be solved for. Substitution of () into () leads (after considerable algebra) to 
„ Q,0, + Q,0„ . , „ Q,0, -Q,0, . 


Q,0 1 +Q 2 0 2 _ Q,© 2 

^>2 , ^2 A i and y ~ r^2 . r^2 


+ fl 2 


ZLL - 1 A 

q 2 +q 2 b 


Q, = D 4 - 6 D 2 co; + col - 2N£)2 + 2N ®o + ® 


Q 2 = 4 D co 0 - ADco 0 - 4N£)<y 0 


The general solution for x > 0 is a linear combination of the complementary and particular solutions in the form 
v(x) = C l e~°” px cos (cox) + C 2 e^ x sin (cox) +C 3 e aipx cos (cox) + C 4 e m/,x sin(<wx) 

+ e~°* nx (Pcos (co 0 x) + Qsin (co 0 x)j 

However, the solution must be bounded as x — > oo so we must have C 3 = C 4 = 0 , leaving 

v(x) = (C, cos(<yx) + C 2 sin (cox)) + e~°’°‘ p ° x (,Pcos(<y 0 x) + (2sin((y 0 x)) (77) 

We can solve for the two unknown constants, C, and C 2 using the boundary conditions at x = 0 provided with 
(62). For this purpose we need several derivatives. These are 

— = -e~ mpx (Cj'cos((Wx) + C 2 sin(®x))-e“' WI (P'cos(<y 0 x) + (2'sin(<» 0 x)) (78) 

dx - / v / 

where 

C; = co(<pC] -C 2 ), C' = co(cpC 2 + C, ) , (79) 

P' = w 0 (cp 0 P - Q) , Q' = a> 0 ((p 0 Q + P) (80) 


as well as 


-jY = +e mpx (C,"cos(®x) + C 2 sm(cox)) + e ° VP " X (,P''cos(<y 0 x) + (2''sin(<» 0 x)) 


C”=(co 2 cp 2 — co) C, -(oo 2 cp + oocp)C 2 , C” = (co 2 cp + co(p)C l +(co 2 cp 2 -co)C 2 
P — ^ co Q (p 0 — co 0 ) P — (^ co 0 (p 0 + cofflQ j Q , Q — ( co 0 (p 0 + co^cp^ ^ P + ( co 0 cp q — co 0 ^ Q 


and finally 


—jY = -e (C”'cos(cox) + C"'sin(cox))-e co °‘ p ° x (P'''cos(o 0 x) + Q'" sm(co 0 x)) 


where 


Cj'"= ^afcp 2 - co)cocp - (a 1 cp + cocj^co\c x -\^a> 2 cp 2 - a)a + ( y a 2 cp + a>(p)a>(p^C 2 (85) 

C 2 " = \ico 2 cp + a>(p)a>(p + (co 2 (p 2 -a>)co^C l -T(® 2 cp + cocp)co -(<y 2 cp 1 -<y)<y^jc 2 (86) 

P — ^(^0^7) — — (^0 *Po ^0^0 ) ^0 J P — ^(^*0 ^0 — ^0 ) ^0 (^0 *Po ^0^0 ) ^0^0 J Q (^ 2 ) 

Q —^(cOqCPq + cDqCp^oOqCP^ +) co f) cpQ — — +m 0 ^? 0 — co 0 ) co^cp^ J Q (88) 

Applying the boundary conditions, we first have zero slope condition at x = 0 , that is 
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( 89 ) 


dv(0) jdx = -C[ - P' = 0 

and the third derivative must be zero meaning there is no jump in the shear load so that 

d i v{Q)/dx i = -C["-P'" = 0 (90) 

These two equations result in C[ = -P ' , C"'= —P”' and thus we obtain the simplification into two equations 

(pcoC i -(oC 2 =~(O 0 {(p Q P-Q) (91) 

and 


[^(<yV 2 - co^cocp -(er cp + -\^m 2 (p 2 - co^co + {a 2 cp + ( 0 (p^acp\c 2 

= -[(® 0 Vo - % - (®oVo + ) ®o ] p ~ [(®oVo - + (®oVo + «Wo ) ] Q 

in two unknowns, Cj and C 2 . Letting 
A n =cpco and A u =-co 

A 2l = \[co 1 (p 2 - co^co(p-[co 2 (p J r cocp^co^ and A 22 = -\^o) 2 (p 2 - co + (w 2 <p + cocp^cocp^ 




and 


B i=~\ 


-[(®oVo -® 0 )®o^o-(®oVo +®o%)®o]' P + [(®oVo -®o)®o+(®oVo + ®o^o)®o^o]2 

and using Cramer’s rule on (91) and (92), the solutions for C 1 and C, are 


C,= 


A 22 B t A l2 B 2 A 2i B\+A n B 2 


A A - A A 
^ 11^22 ^ 21^12 


and C 2 = 


4l^22 ^21^12 

Thus our solution for x > 0 is 

v(x) = e~ mq ’ x (C, cos(®x) + C, sin(®x)) + e”' ! ' >w '’ r (Pcos(® 0 x) + (7sin(fti 0 x)), x > 0 
We also require the bending moments and these depend on the curvature given by 


j2 ~ 

— j- = +e~°” px (Cj"cos(®x) + C 2 'sin(ftix)) + e~ ra " ?, ’’ r (P''cos(cu 0 x) + (7' , sin(tt> 0 x)), x > 0 


where 


C"=[co 2 (p 2 -[afcpA- axp}C 2 , C 2 =[(D 2 (p + a)(p}C x +(a>V 2 -<y)C 2 

P" = ( ®oVo - ®o ) -P - ( ®oVo + ) £? and Q" = {vl<P f) +OJ a (p (i )P + { col <pl -<» a )Q 

We can extend the solution to -oo < x < oo by writing 

v(x) = ^Cj cos(®x) + C 2 sin(®|x|)) + e~“ Wo ^ (Pcos(<» 0 x) + (7sin(® 0 |x|)j, -oo < x < oo 

and 


(92) 

(93) 

(94) 

(95) 

(96) 

(97) 

(98) 

(99) 

( 100 ) 
( 101 ) 

( 102 ) 


j2~ 

= +e~ aq ” i (C"cos(a>x)+ C"sm((o\x\^ + e~ a ° nx (/ > ''cos(fti 0 x) + (7''sin(® 0 |x|)j, -oo < x < oo (103) 

The last two equations represent the key results of the derivation, and allow us to calculate the changes in transverse 
liner displacement, interface pressure and bending moments caused by the frozen in transverse displacement of the 
liner inherited from autofrettage. 

F. Procedure for Calculating Displacements, Moments and Contact Pressures at Zero Vessel Pressure. 

Because of the structure of the fundamental differential equation (62), it is necessary to specify a reference state 
where the liner tension is zero, N, = N = 0 and the frozen-in displacement profile is v 0 (x) . Otherwise, at any 

other value of N, and corresponding value of N , the displacement profile will be v(x) = v(x) + v 0 (x) and will not 
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match the reference state on the right hand side of (62). This would considerably complicate the analysis, since the 
development of the displacement disturbance (depression) from void compaction occurs primarily at some effective 
liner tension, denoted N, c , close to the liner tension at peak proof pressure, denoted A \ . (This is chosen slightly 

different from N, to allow some flexibility in adjusting the reference depression shape, v 0 (x) .) Furthermore, the 
liner yielding causing relaxation of the bending moment and dispersion in the contact pressure around the initial 
very narrow and very stiff region, also occurs at the liner tension, N, , at peak proof pressure. Thus we must 

consider the unloading phase from N, = N, to N, = 0 in order to determine the conditions that exist at peak proof 
pressure compared to the reference state. Once appropriately determined, we can reduce the liner tension from 
A ; = 0 to N, = N, u , which is a negative value corresponding to 0 psi vessel pressure, and estimate the 
displacement, moment and contact pressure conditions that apply there. 

An important point to note is that the measurement of the depression depth profile and maximum value, 8 b u , is 

typically performed at zero vessel pressure using a profilometer or eddy current probe. Thus whatever parameters 
are chosen to determine the reference state and state at proof, they must result in this measured profile. What is 
found is that the liner depression displacements are largest at zero vessel pressure ( A, = A, „ , 8 b = S bu ), only 2 to 

5% smaller at the zero-tension liner state ( N, = 0, 8 b =S b0 ), corresponding to about 2650 psi in a 40 inch vessel, 
and are an additional 12 to 20% smaller at the full proof pressure ( A , = A, , S b = 8 b ). While small, these 
differences are important to consider during modeling for possible buckling. 

Earlier we derived a result for the frozen-in displacement profile at N, = 0 , given by (57). This shape 
corresponded to a sharp line displacement or equivalent force per unit length (transverse to the line of the in-plane 
load) of maximum reference depth magnitude, 8 h 0 , assuming model parameters based on the initially low 

transverse stiffness of the overwrap, E Ol0 , since this is what prevailed when the depression was formed at 
A, = N, . One might question whether the line force ought to be replaced by the dispersed pressure version, but 
this would greatly complicate an already complicated analysis with little benefit since it makes virtually no 
difference to the reference, frozen-in displacement profile, v 0 (x) , at N, = 0 . Also, once the frozen-in deformation 
occurs at A, = A, c , we assume the overwrap stiffness is the same everywhere and now has the stiffened state 
E o,< » E o,u o • Thus, changes to the bending moment profile, when A, is increased or decreased relative to A, = 0 , 

are smooth and involve no discontinuity in slope (as occurs for instance in three-point bending). This calculation 
was the heart of the latest calculation. 

Finally, we must treat the residual bending moment and associated interface pressure disturbance that results 
when the compaction occurs at an effective liner tension A, = A, c » A, at maximum displacement S b 
corresponding to maximum proof pressure. This calculation, performed earlier began with parameters associated 
with the transversely soft overwrap state, E a t 0 , and a sharp line load corresponding to displacement, 8 b . 

Calculating the nominal pressure profile assuming bending in the elastic state, requires obtaining the fourth 
derivative of displacement. In keeping with equilibrium, it was found that the pressure profile involved a force 
component applied as a Dirac delta function along the applied line displacement of magnitude, 8 b . Using a fourth 

order polynomial (46), dominated by a second order parabola, the line-force was redistributed as a parabolic like 
pressure distribution, with the fourth order terms providing continuity in displacement and slope at the boundaries of 
this pressure distribution. Additionally, the nominally elastic bending moment distribution were relaxed to about 
half because of the plastic yielding occurring during autofrettage, which completely blunted half of the stress 
perturbation from bending. Also for higher depression displacements, a plastic limit in bending was assumed to 
occur, restricting the maximum bending moment that can exist at proof. Once the form of the pressure distribution 
was known, the bending moment distribution could be calculated, and a relationship between the maximum 
pressure, width of the dispersed force and maximum bending moment could be calculated. 

Thus the calculation proceeds by first estimating the maximum displacement difference that would occur 
between the reference state A, = 0 and proof state A, = A, , using the solutions developed from (57) to (103). In 

other words, 8 b p -8 b0 -v(o,A /p ) . Once 8 b p is known, we calculate the residual bending moment distribution, 
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M (x) , and interface pressure distribution, p (x) , with reference to this state of maximum liner tension at proof. 
Then using the moment and displacement solution with reference to N, = 0 , we can calculate the change in bending 
moment distribution and pressure distribution in unloading from N, = N, to N, = 0 . Then, as a final step, we can 
reduce the liner tension from N, = 0 to N, = N, u corresponding to zero vessel pressure. The results for the final 
moment distribution, pressure distribution and displacement profiles are 


M u( X ) = M c{ X ) + K P EI (* ; N l,u ) “ ( X; N I.P ) 

(104) 

Pu( x )=p P ( x )+ ^{ v (*;#/,„) 

(105) 

and finally 

v(x) = v 0 (x) + v(x;A,) 

(106) 


The critical conditions for the onset of buckling involve searching for cases with two critical properties: First is 
that the liner pressure goes negative, i.e., the interface has tensile tractions, which cannot be sustained since the 
overwrap is not glued to the liner and the second is that the bending moments have zero points more than two inches 
apart spanning the negative pressure region. When both these conditions exist, it is possible for buckling to occur. 

V, Results of Model Calculations on Potential for OMS Liner Buckling 

We now present graphical results from the above analysis for a 40 inch Kevlar/epoxy spherical COPV, with titanium 
liner of thickness, t, = 0.106 in. , Young’s modulus, E, = 16.5x10'’ psi , and Poisson’s ratio, v t = 0.342. The 
maximum possible value for N \ during autofrettage is determined from the liner yield stress, a, = 128,000 psi as 
N, p = <7 l y t[ = 1 3,570 lb/in. The liner compressive load at zero pressure after autofrettage is nominally 
N lu =-9,172 lb/in. In the mathematical calculations, a reference state of zero liner tension, N,=0 lb/in., was 

chosen, which happens to corresponds to 2,650 psi pressure in the vessel after autofrettage. Choice of this reference 
state greatly simplified the mathematical analysis, but has no effect on the predicted results and conclusions. 

Based on the analytical results for depression displacement, v (x) , at liner tension, N, , after autofrettage, Figure 
5 is the predicted profile from a 40 mil depression in 40 in. liner corresponding to the reference state at zero liner 
tension, N,=0 lb/in. , at a vessel pressure of 2,650 psi following autofrettage. The depression depth is measured 
relative to an equilibrium radial deformation state for a perfect sphere, balancing the vessel pressure and the liner 
and overwrap tension. In the formation of this depression during initial pressurization, the through-thickness 
stiffness of the overwrap was taken as E 0 0 = 25,000 psi , representing the pre -proof state with voids (and 
consistent with the result from the special OMS vessel tested in 1977 at JSC). 

For the unloading phase from the peak proof pressure of 6,520 psi, and peak liner tension N, = 13, 570 lb/in. the 

through thickness stiffness of the overwrap was increased to E ot =250, 000 psi, due to overwrap compaction 

consistent with experiments at WSTF using eddy current probes. The resulting depression width (half-wavelength) 
of about 5 inches matches well the experimental width measured by laser profilometry in Figure 3. Note that under 
the full tension of the proof pressure, the depression depth of about 37.75 mils is less than the 40 mil reference depth 
at zero tension, and under full liner compression at 0 psi vessel pressure, the peak displacement is 42.25 mils, 
slightly larger than the 40 mil reference depth. These differences are relatively small because of the greatly 
increased through-thickness overwrap stiffness once void compaction is complete. 

Note also that the incremental profiles representing changes in displacement measured relative to the reference 
state have much shorter half- wavelengths than the originally formed depression, again because of this greatly 
increased overwrap stiffness. This mismatch in half-wavelength turns out to be a very important feature in 
controlling the onset of buckling. 
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Figure 5. Predicted characteristics of the displacement profile in an OMS liner due to a nominal 40 mil 
depression and shown at the reference state and the states of proof pressure and zero pressure after proof. 
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Figure 6. Predicted characteristics of the interface pressure profile in an OMS liner due to a nominal 40 mil 
depression and calculated at the proof pressure and at zero pressure after proof. 

Figure 6 shows the predicted interface pressure profile in an OMS vessel from a 40 mil (0.040 in.) depression 
corresponding to the displacement profiles in Figure 5. This includes the nominal or average interface pressure bias, 
p 0 =956 psi at zero vessel pressure, which results from radial interference of the overwrap and liner induced by 
plastic yielding of the latter during autofrettage. Also shown is the state of the interface pressure profile at the peak 
proof pressure of 6520 psi, and N lp = 13,570 lb/in. as well as the changes in interface pressure from the proof 
pressure to the reference state of zero liner tension (at 2650 psi), and from this reference state down to zero vessel 
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pressure. Note that the delta function spike in pressure due to the idealized line load at proof pressure evolves into a 
more dispersed, approximately parabolic profile from |x| < 0.65 in. A key observation is that the interface pressure 

drops to about 400 psi over close to a 1.3 inch region centered over the depression. Still there remains a reserve of 
interface pressure to maintain firm contact between the overwrap liner thus maintaining elastic stability. 

Figure 7 shows the calculated bending moments associated with the 40 mil depression under the conditions of 
Figures 5 and 6. Shown are the bending moment distribution at the peak autofrettage pressure as well as the changes 
in bending moment upon reducing pressure to the reference state of zero tension in the liner at 2650 psi, and then 
from this reference state down to the fully unloaded state of zero vessel pressure. We sum the two latter 
components to obtain the total change in bending moment distribution from proof to zero vessel pressure, as well as 
the final bending moment distribution upon summing all components. 

An important observation in Figure 7 is that magnitude of the change in maximum bending moment coming 
down from the peak proof pressure differs very little from the final magnitude. Thus, unlike the situation with the 
interface pressure, the magnitude of the residual bending moment at proof has little influence on the final bending 
moment at zero vessel pressure. The zero moment hinge points are considerably wider than 2.0 inches meaning that 
significant benefit of overwrap support exists outside of rage the reduced interface pressure. Flowever the final 
bending moment is large enough to impinge on the limitations imposed by the Bauschinger effect, and even if the 
final bending moment might increase with an increase in depression depth, say to 50 mils, softening in the liner is 
likely in bending as reflected by a reduced secant modulus seen in Figure 4. 
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Figure 7. Predicted characteristics of the bending moment profiles in an OMS liner due to a nominal 40 mil 
depression and calculated at the proof pressure and at zero pressure after proof. 

Figure 8 shows the change in interface pressure that results at 40 mils depression depth when the effective liner 
stiffness (secant) is reduced 15% from 16.5 Msi to 14.0 Msi. This change is enough to reduce the interface pressure 
to zero over a substantial region spanning the depression, thus resulting in critical conditions for liner separation. 
The same loss in secant stiffness of the liner material in bending leads to further increases in elastically calculated 
bending moment, as shown in Figure 9. This signals an unstable state where yielding due to the Baushinger effect 
on depressurization does not lead to bending moment relief but rather runaway deformation signaling liner collapse. 
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Figure 8. Effect on interface pressure of reducing the liner stiffness 15% from 16.5 Msi to 14.0 Msi, when the 
nominal depression depth is 40 mils. 
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Figure 9. Effect on the bending moment of reducing the liner stiffness 15% from 16.5 Msi to 14.0 Msi, when 
the nominal depression depth is 40 mils. 

The Bauschinger limit imposed in Figure 7 may be conservative, so we consider the sensitivity of these results to 
changes in the depression depth. Figures 10 and 11 show the effects of changing the depression depth for the cases 
in Figures 5 through 7 and returning to a liner stiffness of 16.5 Msi. Clearly it is highly unlikely for the liner to be 
stable at depression depths above about 55 mils irrespective of the details of Bauschinger yielding effects. 
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Figure 10. Effect on the nominal interface pressure distribution in an unpressurized vessel of changing the 
nominal depression depth over a range from 20 mils to 80 mils. 
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Figure 11. Effect on the nominal (elastic) bending moment distribution in an unpressurized vessel of 
changing the nominal depression depth from 20 mils to 80 mils. 
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VI. Concluding Comments 

A key observation is that the magnitude of the change in maximum bending moment when depressurizing from 
the peak proof pressure differs little from the final magnitude. Thus, unlike the situation with the interface pressure, 
the magnitude of the residual bending moment at proof has little influence on the final bending moment at zero 
vessel pressure. The zero moment hinge points are considerably wider than 2.0 inches meaning that significant 
benefit of overwrap support exists outside of range of the reduced interface pressure. However, the final bending 
moment is large enough to impinge on the limitations imposed by the Bauschinger effect, and even if the final 
bending moment were to increase with increased depression depth, softening in the liner is likely in bending as 
reflected by a reduced secant modulus. 

The main conclusion is that the liner is likely to stable and not buckle from depression depths of up to about 40 
mils. Fortunately this depth is double the depression depths of about 20 mils that have been observed using laser 
profilometry, and are believed a characteristic of the winding process. Figures 10 and 11 indicate that a nominal 
depression depth of 20 mils (actual about 21 mils) presents no hazard. Above 40 mils, the Bauschinger effect 
becomes extremely important and destabilization and buckling becomes increasingly likely depending on the details 
of depression formation during the proof pressurization. Whatever the mechanical property variations for a 
particular vessel, it appears almost certain that destabilization and buckling will occur for depression depths after 
proof beyond 55 mils, and beyond 40 mils there is considerable risk. 
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